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ABSTRACT 

The findings of a nine orbit calibration plan carried out during HST Cycle 15, to fully determine the 
NICMOS camera 2 (2.0/zm) polarization calibration to high accuracy, are reported. Recently Ueta et 
al. and Batcheldor et al. have suggested that NICMOS possesses a residual instrumental polarization 
at a level of 1.2 — 1.5%. This would completely inhibit the data reduction in a number of GO pro- 
grams, and hamper the ability of the instrument to perform high accuracy polarimetry. We obtained 
polarimetric calibration observations of three polarimetric standards at three spacecraft roll angles 
separated by ~ 60°. Combined with archival data, these observations were used to characterize the 
residual instrumental polarization in order for NICMOS to reach its full potential of accurate imaging 
polarimetry at p ks 1%. Using these data, we place an 0.6% upper limit on the instrumental polar- 
ization and calculate values of the parallel transmission coefficients that reproduce the ground-based 
results for the polarimetric standards. The uncertainties associated with the parallel transmission 
coefficients, a result of the photometric repeatability of the observations, are seen to dominate the 
accuracy of p and 9. However, the updated coefficients do allow imaging polarimetry of targets with 
p « 1.0% at an accuracy of ±0.6% and ±15°. This work enables a new caliber of science with HST. 
Subject headings: Astronomical Instrumentation. Data Analysis and Techniques. 



1. INTRODUCTION 

Polarimetry is a powerful observational tool that aug- 
ments and complements the capabilities of imaging, pho- 
tometry and spectroscopy. While the latter allow the 
determination of spatial distributions, chemical compo- 
sition and dynamics, polarimetry lets us observe the na- 
ture of magnetic fields, object orientation, scattering, 
and the properties of interstellar particles in general. It 
also allows us to probe the nature of emission mecha- 
nisms (e.g., synchrotron/thermal), and to investigate the 
geometry of unresolved sources. 

The value of polarimetry to elucidate the physical 
nature of light-scattering particles and their environ- 
ments in astrophysical systems has been repeatedly 
demonstrated in, for example, active galac tic nuclei 
(jAntonucci fc Miller! Ii~985l : ICapetti et allll995[). gal actic 
and extra-ga lactic magnetic fields (jWi elcbinski 199j|), ra- 
dio galaxies (iFeretti et al.lfl9 98: Vcr net et al.l [2001). and 
circumstellar disks (Ku hn et al.l 2001). However, because 
of "beam depolarization" , results derived from polari- 
metric imaging observations may be biased when ob- 
tained at low spatial resolution. For example, compact 
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(but spatially resolved) objects that exhibit circular sym- 
metry or bi-conical structures, like active galactic nuclei 
(AGN) and stars surrounded by orbiting circumstellar 
dust, contain polarized elements that partially reduce, 
if not fully cancel, each other when spatially averaged. 
If the object is unresolved, then the polarization vectors 
sum to zero; in any case, beam depolarization will always 
underestimate the intrinsic polarization. 

NICMOS is the only near infrared (NIR) instrument 
capable of the high resolution, high fidelity polarime- 
try needed to examine the scattering geometry and ma- 
terials, in detail, for many types of astronomical ob- 
jects. To date, there have been a multitude of stud- 
ies carried out using NICMOS d irect imaging polarime- 
try (e.g. ICapetti et all l2000t iTadhunter etall 120001; 
I Silber et all 120001: [Schmidt et all 120021: I Simpson et~al 
20021 ISu et al.l l2003t iMeakin et all 120051: ISimpson et al 
|2006|) and there are now more ambitious programs to 
measure the characteristics of AGN and regions within 
circumstellar protoplanetary debris disks, that typi- 
cally exhibit polarizations of less than 5% (ISmith et al.l 
l200l iSchneider fc HST/GO 10177 Team! I2005h . How- 
ever, there have been recent reports of a systematic 
residual instrumental p olarization of ~ 1.5% in some 
objects (|Ueta et al.l 120051 ) . which although acceptable 
for highly polarized targets, seriously compromises stud- 
ies of low polarization objects. A comprehensive study 
of this residual polariz ation has been carried out by 
iBatcheldor et al.l ((2006, hereafter B06). Using all of the 
available standard polarization data in the archives, B06 
indeed find a residual excess polarization of w 1.2% in 
the NIC2 camera. The extent of this excess remains un- 
known for NIC1. 

Unfortunately, prior to the new calibration reported 
in this paper, since the adven t of the NICMOS cooling 
system, NCS (|Schneiderll2004r i. installed during SM3B in 
2002, observations of only one polarized and one unpolar- 
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TABLE 1 
Details of Observations 



Target 


DATE-OBS 


ORIENTAT 


Camera 


T (s) 


VR 84c 


2007-05-03 


182.2° 


NIC1 


3.6 




2007-08-03 


299.3° 








2006-12-15 


59.2° 








2007-05-03 


181.4° 


NIC2 


14.0 




2007-08-03 


on a O 

298.6 








2006-12-15 


58.4° 






VSS VIII-13 


2006-08-04 


180.3° 


NIC1 


1.6 




2006-06-25 


114.3° 








2007-04-11 


54.3° 








2006-08-04 


179.6° 


NIC2 


7.0 




2006-06-25 


113.6° 








2007-04-11 


53.6° 






HD331891 


2006-08-27 


273.8° 


NIC1 


1.6 




2006-07-01 


338.4° 








2007-04-15 


38.3° 








2006-08-27 


273.1° 


NIC2 


7.0 




2006-07-01 


337.6° 








2007-04-15 


37.5° 







Note. — Each observation was performed at three different 
roll angles and for the POLO, POL120 and POL240 polariz- 
ers. In total, NIC2 was dithered around 36 pointings per orbit, 
with 12 pointings per polarizer. The DATE-OBS parameter 
gives the observation date in YYYY-MM-DD format. "ORI- 
ENTAT" refers to the celestial position angle of the image +y 
axis (degrees east of north). "T" gives the exposure times in 
seconds. 



ized standard have been obtained, and each at only two 
celestial orientation angles. Hence, it has been impossi- 
ble, with those limiting data, to entirely characterize the 
NICMOS residual polarization. Polarization measure- 
ments at three well separated roll angles are required to 
remove the dependence of the measured Stokes parame- 
ters on the relative transmission of each of the filters. 

This paper reports the findings of a nine orbit HST 
calibration plan (Cycle 15), to fully investigate the 
residual polarization in NICMOS. NICMOS carries two 
cameras with polarimetry optics designated NIC1 and 
NIC2 (a third camera for wide-field and grism imag- 
ing has no polarimetry capability). NIC1 provides an 
ll"xll" field-of-view (FOV) with an image scale of 
43.1 milli-arcsecs (mas) per pixel and 1.045/im broad- 
band (0.475Atm FWHM; R = AA/A = 45%) po- 
larizing filters. NIC2 (the primary focus of this pa- 
per) provides a 19"3x 19'.'2 FOV with 75.8 mas pix- 
els and a medium bandwidth polarimetric passband 
of 1.994/im (0.202^m FWHM; R = 10%). Both 
optical channels critically sample the respective point 
spread functions (PSFs) in their polarimetric passband. 
For additional d e tails of the NICMOS instrument see 
iThompson et al.l 111 9981) and the NICMOS Instrument 



Handbook ( Barker et al.ll2007l ). The two cameras have 
non common path optics from the instrument's field di- 
vider mirror assembly to their respective re-imaged focal 
plane detectors. As polarimetric imaging in the two cam- 
eras is (uniquely) carried out at different wavelengths, 
there is no reason to expect the same level of residual 
instrumental polarization in each camera. With this in 
mind we carried out observing plans to push the NIC2 
calibration to its instrumental limit, given the intrin- 
sically superior pol arimetric perform ance of the NIC2 
camera over NIC1 (|Hines et al.l |2000| ). Due to the na- 
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Fig. 1. — An example from one orbit of the NIC2 pointings 
used to determine the effects of the IPRFs, quadrant dependence 
and persistence. POL0L, POL120L and POL240L observations are 
in the inner, middle and outer "diamonds" of points respectively. 
The "A", "B" and "C" positions are used to index the dithered 
pointings of each polarizer. The POL0L, POL120L, POL240 and 
A, B, C markings, that are both applicable to each quadrant, are 
used to define the pointing naming convention. For example, the 
lowest pointing on the left of the array would be POL240L-Q3-C. 

ture of NIC1, we confine the results from this camera 
to Appendix [A] In § [2] we describe the observing strat- 
egy employed and in § [3] we explain the data reduction 
procedures. A detailed examination of the photometry 
is completed in § 13.11 The results and methods are pre- 
sented in §2] before being discussed in §[3] and concluded 
in § [6j A recommended observing plan for NICMOS 



imaging polarimetry with NIC2 is included in § 15.11 

2. OBSERVATIONS 

The NIC2 observations were designed to be totally 
comprehensive in the sense that they were not only made 
to address the NICMOS residual polarization, but sev- 
eral other key issues including possible dependence with 
detector quadrants, source color dependence, the inter- 
pixel response functions (IPRFs) and latent image per- 
sistence. The details of the observations are summarized 
in Table □ 

Image persistence (or latent image decay) can con- 
tribute significantly to photometric errors in sequential 
exposures of bright targets on the same pixel locations 
on the focal plane array. The same is also true for cos- 
mic ray hits, therefore all observations were scheduled 
in orbits that were not impacted by the South Atlantic 
Anomaly. As NICMOS polarimetric analysis derives in- 
tensity differences from measures obtained in each of its 
three camera-specific polarizing filters, systematic dif- 
ferences due to image persistence, or other instrumen- 
tal causes, should be circumvented or minimized to the 
greatest degree possible. Persistence can be overcome 
with an effective image dithering strategy. The diffrac- 
tion spikes arising from the HST secondary mirror sup- 
port in the J25T+NICMOS PSF can contain a large num- 
ber of photons, so a precise dithering pattern has to be 
carefully chosen. In our calibration observations, custom 
dither patterns were created for each camera. 

In the absence of super-dark calibration frames, 
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TABLE 2 
Details of Polarimetric Standards 



Target Type m v rrij mj 0.55/jm 1.05/jm 1.21^tm 2.00/xm 2.04/im 
(p%,0) 

VR 84c (WKK F7) B5V 10.3 9.34 9.13 5.98,118° 4.04,118°* 3.19,118° 1.25,126°* 1.19,126° 
VSS VIII-13 K1III 12.7 9.65 8.68 2.25,102° 2.49,99°* 2.21,96° 0.91,102°* 0.86,102° 
HP 331891 A4III 9.3 8.80 8.72 0.04,n/a — -j_ 

Note. — Basic data for the observed polarimetric standards. Two polarized and one un-polarized star 
(HD331891). Values marked with "*" are those that have been corrected using the Serkowski curve. There is 
no method to precisely determine 6 at 1.05/^m and 2.00/xm so the average values between 0.55/xm and 1.21^tm, 
and the values at 2.04/mi are used respectively. More details for each target are given in the text. 



that replicate to very high fidelity the behavior of 
the NICMOS detectors as specifically clocked for any 
given science observation, both quadrant dependent and 
global bias offsets often arise in basic (pipeline) cal- 
ibrations (both are frequently seen when calibration 
database/library or synthetic dark frames are employed). 
However, during this calibration we acquired contempo- 
raneous dark frames to mitigate the possibility of residual 
systematics in the multi-accum derived count rate im- 
ages. These darks were matched to the detector clocking 
modes used for our targeted observations, and executed 
during the occulted period of each orbit. None-the-less, 
to assess the possible impact of "pedestal-like" effects of 
other instrumental origins, the same observing sequence 
was repeated with the target repositioned in each of the 
four detector quadrants. Flat-field artifacts can be en- 
hanced (rather than reduced) in amplification with im- 
properly biased (subtracted) dark reference frames. 

The observations in each quadrant were dithered in 
order to mitigate the effects of the IPRFs and to sam- 
ple around defective pixels. The IPRFs can alter the 
count rates from a target simply through different align- 
ments (centerings) with respects to the pixel center at 
each pointing. The pointing offsets cover a phase spac- 
ing of -1/3, and +1/3 pixels in order to best tile the 
IPRFs for the dither pattern employed. In addition, the 
dither step sizes are greater than the radius of the 3 rd 
Airy PSF minima (0'.'58, 7.6 pixels). The orientation of 
the dither pattern was rotated, with respect to the detec- 
tor rows and columns, by 22.5° to remove the possibility 
of overlapping persistence from the PSF diffraction spikes 
of previous pointings. Figure [T] shows an overlay of the 
NIC2 pointings from within one orbit. 

These observations were repeated for each of the three 
chosen polarimetric standard stars and at three sepa- 
rate spacecraft roll angles, i.e., three celestial orienta- 
tion angles of the camera aperture on the sky. This 
enables a unique determination of any residual instru- 
mental polarization, above the detection floor, relative 
to the equatorial reference frame. More crucially, this 
multi-orientation technique means that the polarimet- 
ric analysis derived from the three polarizers at each of 
the field orientations can be completely and uniquely de- 
coupled and checked for consistency. This method also 
removes the dependence of the measured Stokes parame- 
ters on the relative transmission of each of the polarizing 
filters and, at the same time, any dependence of the in- 
strumental polarization on the detector's sensitivity to 
polarized light. 

Light scattering within the instrument could, at least 



in part, be largely responsible for the observed residual 
instrumental polarization, but light scattering is wave- 
length dependent, so standard stars with intrinsically 
different spectral energy distributions were chosen. The 
details of the three polarimetric standards picked for this 
calibration are summarized in Table [2] As HST was op- 
erating in two gyro during HST Cycle 15 (the epoch of all 
the calibration observations - see Table [1]), the dominant 
factor in the selection of this sample was the availability 
of the three separate roll angles required to allow po- 
larization measurements to be derived from the single 
polarizers. 

2.1. Calibration Targets 

All targets were sufficiently bright so that the PSF cen- 
tral pixel approached the full-well depth near the end of 
the non-destructive sequence of readouts in each multi- 
accum exposure. This maximized the photon signal to 
noise ratio with the greatest observing efficiency. In ad- 
dition, as the exposure times were short (Tabic Q}, the 
target brightnesses ensured each could be observed mul- 
tiple times within its dither pattern during one orbit. 

VR 84c (WKK F 7, Cha DC F7) was chosen as it has 
a significant percentage polarization in both the J and 
K bands, and as it had been observed in the previous 
pola rimetric calibrations o f NICMOS; it serves as a con- 
trol. IWhittet et all (|1992j) give a K-band polarization of 
1.19 ±0.01% at 126 ±4°. 

The second polarimetric standard VSS VIII-13 (R CrA 
DC No. 13) was also chosen for its percentage polariza- 
tion. However, it is also found in a different associ- 
ation from VR 84c making it more suitable as a cal- 
ibration target. Within any given association, mea- 
sured stellar polarization angles tend to be similar. This 
implie s a polarization pro cess internal to the associa- 
tion. Whittet et al.l (|1992t ) give a K-band polarization 
of 0.9 ±0.2% at 102 ±5°. 

The final standard, HD 331891, was chosen as it was 
found to have an insignifi cant amount of polariz ation at 
0.55^m (0.04 ±0.02%) bv lTurnshek et all (|1990l ). It too 
was observed in previous calibration projects. There are 
no data of this target at 2.00/im(previous to this study), 
therefore we have to cautiously assume that it remains 
unpolarized at this wavelength. 

Comparison of NICMOS and ground-based polariza- 
tion measurements must be band-transformed through 
the known wavelength dependen ce of interstellar po- 
larization (Serkowski e t al.l Il975h . The ground-based 
results are derived from observations at 1.21/xm and 
2.04/im and are corrected to 1.05/im and 2.0/im using the 
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Fig. 2. — For each target, the encircled count rate summed over all dither positions, for each polarizer and target field celestial orientation, 
is shown. The inflexions due to the PSF are clearly visible, as are some deviant profiles (e.g., HD 331891 at 37.5°). In all cases 95% of the 
total flux is with a radius of 20 pixels, 99.9% is at 50 pixels. 



"Serkowski curve" (|Serkowskilll97llWilking et al.lfl98l ). 
These values are listed in Table [5] and are used in this 
study as the ground-based comparison values. We also 
propagate the uncertainties of each ground-based mea- 
surement in our analysis. 

3. DATA REDUCTION 

The "raw" multi-accum exposures were used to cre- 
ate individual instrumcntally calibrated count rate im- 
ages for each target pointing. Calibration steps in- 
cluded: "DC bias" offset corrections with respect to 
the contemporaneously obtained reference dark frames, 
"super-zero" signal accumulation compensation, read- 
by-read dark subtraction, read-by-read linearity correc- 
tion, comic-ray identification and reparation, saturation 
threshold determination, flat-field correction,weighted- 
least squares count rate solution, bad pixel reparation 
by 2D interpolation of good neighbors and the removal 
of the telescope (+sky) thermal background. Before con- 
tinuing, each image was checked for possible latent image 
persistence resulting from the previous image. No image 
latency was found in any of the individual multi-accum 
exposures. 

3.1. Photometric Analysis 

It is essential to individually evaluate the fidelity of 
each pointed observation. An ideal imaging detector 
would produce precisely the same photometry with the 
target located at any arbitrary position on the grid of 
detector pixels, and would make polarimetric analysis 
fairly straightforward. Polarimetry, by its very nature, 
is highly sensitive to small variations in the photometry, 
and in any instrumental system there are many factors 
that can compromise the data. Sources of small varia- 
tions in the photometry, that do not directly result from 
the intrinsic polarization of the target, and that have 
not already been accounted for, must be identified and 
removed. 

Aperture photometry was performed at each pointing 
using the digiphot package within IRAF 8 . Radially in- 
cremented circular apertures were placed on the target. 

8 IRAF is distributed by the National Optical Astronomy Ob- 



The (sub-pixel) centers were determined using the inten- 
sity weighted means of the profiles in the x and y di- 
rections. The encircled count rates through the radially 
increasing apertures were then extracted. 

The summed extracted NIC2 radial profiles from the 
same polarizer at the separate position angles are pre- 
sented in Figure [2j From this a number of effects are ob- 
served. Firstly, the Airy pattern component of the PSF 
is seen in the radial profile inflections. Secondly, there 
are some profiles that are inconsistent with the majority 
due to non-repeatability affects in the photometry be- 
yond those mitigated by our observing strategy and im- 
age calibration. Finally, there is a noticeable spread in 
the count rates; more than one expects from polarization 
alone. 

Each individual radial profile was compared to the av- 
erage profile from the 12 individual dithers. The offset of 
each individual radial profile from the average profile was 
determined as a function of radius. Any encircled radial 
profile that showed more than a 2a deviation from the 
average profile was flagged and inspected. Figure [3] spa- 
tially plots the positions of the flagged profiles using the 
filled squares. In all, 10% of the profiles show deviations 
of > 2cr from the original average profile. 

About 40% of the deviant profiles (4% of all the point- 
ings) are the result of the individual target pointings 
placing the center of the stellar PSF on uncorrectable bad 
pixels (excessive grot, i.e., pixels with reduced through- 
put due to particulate contamination on the detector 9 ). 
In this case it is impossible to recover the lost flux and 
these pointings are discarded from further study. These 
profiles are typically clustered around the Ql-A posi- 
tions. 

As expected, in ~ 20% of the deviant cases the point- 
ings have fallen close to the coronagraphic hole. The 
profiles are clustered around the Q2-A positions. With 
time, the precise position of the hole wanders by about 
one pixel. However, its position (x,y) is typically 73.5, 
213.5 with a soft edge that can extend to a radius of 7 

servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 

the National Science Foundation. 

9 http: / / www.stsci . edu /hst / nicmos / documents /isrs /isr_99_008 . pdf 
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Fig. 3. — Identifying NIC2 deviant profiles. The open cir- 
cle (r=0'.'3, 4 pixels) marks the position and geometrical size of 
the coronagraphic hole projected onto the detector. However, the 
afocal imprint of the hole can affect 2.0/xm circum-hole photome- 
try (at our 2cr rejection level) to a distance of ~ 0'.'6. "+" signs 
show the positions of all the pointings from all 9 orbits. The filled 
squares arc profiles that show > 2<r deviations at all radii. The 
open squares are profiles that only show > 2a deviations within a 
radius of two pixels. 

pixels. Data that were affected by the presence of the 
coronagraph were excluded. 

The remaining flagged profiles show > 2a deviations 
at radii of < 2 pixels. As the NIC2 polarizers have an ef- 
fective wavelength of 1.995^tm, and the effective primary 
aperture is 2.281 m (as defined by the NIC2 pupil mask), 
the theoretical PSF core diameter is 0'.'22 or 2.91 pixels 
(PSF FWHM = 0"185, 2.44 pixels); NIC2 critically sam- 
ples the PSF. As these profiles are only deviant within 
the diffraction limit of the instrument, and are within 2a 
of the average profiles at radii > 2 pixels, i.e., the effects 
are negated by increasing the aperture size, they are not 
excluded from further study. 

The processes of flagging and inspecting deviant pro- 
files was repeated until all profiles fell within 2a of the 
average. At each iteration, excluded profiles were re- 
moved from the original average. In the worst case, only 
three profiles (from 12) were rejected. 

3.2. Achievable Accuracies 

In polarimetric analyses, the accuracies that are 
achievable are determined by the signal-to-noise ratio 
(S/N) of the data. Therefore, in order to estimate the 
effects of background noise, the twelve NIC2 images in 
each polarizer were combined, after precise astrometric 
registration, into a single median image. In all cases, im- 
ages were first aligned using the spacecraft pointing in- 
formation provided in the .fits headers of the raw files. 
This co-registered the images with a relative precision of 
a few tenths of a pixel. Subsequent "fine" alignments of 
the inter-visit images were performed using an apodized 
bi-cubic re-sampling of the stellar PSF cores onto a 32 x 
(in both x and y) larger grid. The re-sampled individ- 
ual images were then shifted to a common fiducial posi- 
tion by the difference in centroids found iteratively with 
least-squares 2D Gaussian profile fitting. This process 
was repeated three times and converged all re-registered 
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Fig. 4. — Achievable polarimetric accuracies as a function of 
(S/N) j. and intrinsic target polarization (p(%)). The solid colored 
curves show the increasing uncertainties in p, given by the inset, as 
p approaches zero for a given (S/N) fc . The dashed horizontal lines, 
from top to bottom, are the peak per pixel (S/N)j- for VR 84c, 
VSS VIII-13, and HD331891 respectively, averaged across all three 
polarizers. For example, using the data collected on VR84c, we 
expect to reach accuracies of p = 1.0±0.2% or 0.5±0.4% depending 
on the actual intrinsic polarization of VR84c. Polarizations below 
0.35% will be consistent with zero as the uncertainties are ±0.6%. 

image centroids to a few thousandths of a pixel. 

The pixel-to-pixel noise in the combined polarimetric 
total intensity images (Stokes /) was found by measur- 
ing (in each image) the standard deviation in the "back- 
ground" region. The region was close enough to the tar- 
get that all 12 backgrounds from the individual images 
contributed to the combined image, but sufficiently far 
out so no flux from the target itself contributed to the 
background. A 151 x 151 pixel sub-array region, cen- 
tered on the target was used. Flux from the diffraction 
spikes was masked. To find the region beyond which 
the target did not contribute to the background, square 
masks, increasing in size iteratively, were placed in the 
photometric aperture. As the mask was incremented (in 
2 pixel lengths and widths) from 90 to 100 pixels, ear- 
lier (smaller mask size) standard deviations ceased to de- 
crease, and both the median and mean were consistent 
with zero. Therefore, a square mask of 100 x 100 pixels, 
plus a diffraction spike mask, were used and noise statis- 
tics were computed on 10660 pixels. In all cases, the 
peak pixels (with instrument intensities of several thou- 
sand counts per second) typically have a la noise of 0.1 
to 0.2 counts per second; the per pixel data is photon 
noise dominated; 

As shown bv I Sparks fc Axon! (|1999L hereafter SA99), 
the S/N ratio averaged across the three polarizers 
((S/N)fc) times the required polarization degree (p) de- 
termines the accuracy to which polarimetry can be per- 
formed (ap) in the photon noise dominated regime. 
EquationO] taken from SA99 (their section 7.4), demon- 
strates the relationship between (S/N)^ and a p . It does 
not take into account any instrumental polarization, nor 
any differences in the PSFs between polarizers, and as- 
sumes that the polarizers are perfect. 



Iog 10 (op/p) = -0.102 - 0.9898 Iog 10 (p(S/N) fc ) (1) 

Given p and S /N we can use Equation Q] to calculate 
the theoretical achievable polarimetric accuracies in the 
absence of non photon noise dominated statistics, as il- 
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Fig. 5. — The radial variation of S/N CT in each target. S/No- 
is defined by the values of the average profiles divided by the ltr 
dispersions derived from the individual profiles (see text for more 
details). The three thick vertical dashed lines show the theoretical 
radii of the 1 st , 2 nd and 3 rd Airy minima arising from the diffraction 
pattern of the centrally obscured primary mirror. 

lustrated in Figure [4j It can be seen (by looking at the 
intercepts of the highest dashed horizontal line and the 
±0.2% curve) that the observations theoretically allow 
measurements of p = 1.0 ± 0.2% for VR84c. 

However, as a single NIC2 pixel spatially under- 
samplcs the PSF through the polarizers by 59%, accurate 
polarimetry cannot be done on a single pixel; the image 
plane of NIC2 better than critically (Nyquist) samples 
the FWHM of the PSF. In addition, there are temporal 
instabilities in the PSF from spacecraft breathing and 
pupil mask shifts (the three polarizers are used non si- 
multaneously). For point-source polarimetry, some spa- 
tial binning beyond an instrumental resolution element 
must be done to improve the accuracy of the polarime- 
try otherwise affected by the aforementioned instabilities 
and measurement limitations. 

So far we have assumed that the intrinsic pixel cal- 
ibration uncertainties are largely mitigated through an 
effective image dithering strategy, as was implemented 
for NIC2 in this calibration program. However, there re- 
mains a clear dispersion in the photometry after dither 
combination. As expected, this dispersion varies with 
radius because the diffracted (and scattered) energy ex- 
tends very far from the central pixel (Figure [2|). In addi- 
tion, for any fixed radius, the dispersion is also separately 
affected by target position on the array due to instrumen- 
tal affects (e.g., imperfect flat fielding and IPRFs that, 
in detail, differ among detector photodiodes/pixcls). 

In the point source case we use the encircled energy 
curves, from the incrementing radius aperture photom- 
etry, to follow the variation in the S/N (and therefore 
the achievable accuracy). Small temporal variations in 
the PSF structure of point sources will be largely miti- 
gated by polarimetric analyses done with measurements 
on sufficiently large target "enclosing" photometric aper- 
tures. The optimal size of the photometric measurement 
aperture will be large enough to minimize the aggregate 
affects of changes at the pixel level, but not so large as 
to include data in the noise regime that is not photon 
dominated. Therefore, the S/N from dithered images is 
defined as the ratio between the mean of the individ- 
ual sums of per-pixel measures in each image, and the 
la (standard deviation) of the sums about that mean 



Polarizer 


<f>k (°) 


Vk 


h 


t k 


POL0L 


8.84 


0.7313 


0.1552 


0.8779 


POL120L 


131.42 


0.6288 


0.2279 


0.8379 


POL240L 


248.18 


0.8738 


0.0673 


0.9667 



Note. — The previous polarization coeffi- 
cients as derived from the Cycle 11 calibration 
program. 



(S/No-). With this definition we can see the polarimet- 
ric accuracy possible considering instrumental and tele- 
scope systematic effects (quadrant dependence, IPRF, 
PSF, telescope breathing, etc). 

Figure [5] shows the NIC2 S/N^ as determined from the 
lcr dispersions about the mean in the incremental-radius 
aperture photometry. The typical S/N CT that encircles 
the 2 nd Airy maxima is 300. From Figure [J] we see that 
a polarization measurement of p = 1.3 ± 0.2% is achiev- 
able. This is consistent with the pixel-to-pixel (S/N)*. 
accuracy. The S/N CT that encircles the 1 st Airy maxima 
gives a possible polarization measure of p w 1.2 ± 0.3%, 
inside the first Airy minimum we seep ss 1 . 1 ± 0.5% is at- 
tainable. Apertures with radii of 7.6 pixels are optimum. 
Inside this radius S/N CT is dominated by changes at the 
pixel level (large variations with radius), outside this ra- 
dius photon noise begins to dominate (S/N ff smoothes 
out and turns over). At even larger radii, the target flux 
will strongly decline so that instrumental noise will dom- 
inate. 

The achievable accuracies are optimized by using a 
photometric aperture size of 7.6 pixels (Of! 58). More- 
over, in § 14.11 we will demonstrate that by re-calibration 
of the polarizer transmission coefficients, using our new 
data set, we can bring the observed values of p and the 
celestial orientation of polarization (9) in agreement with 
the ground-based measurement of the polarized standard 
stars. However, as demonstrated by Figure [H in the case 
of the un-polarized standard, we cannot expect to deter- 
mine p = 0%, but merely a polarization that is consistent 
with zero from within the errors. 

4. POLARIMETRIC ANALYSIS 

The determination of p and 9 ha ve been specificall y 
addressed for NIC2 previously (B06; iHines eTaTI [2000l) . 
The coefficients presented in Table |3] (the best available 
previous determination without the addition of our newly 
acquired calibration data) are those determined from 
the Cycle 11 NICMOS polarimetric calibration program 
(|Hinesll2~002f ). We briefly revisit the linear technique for 
the case of three non-ideal polarizers below; the case of 
three "ideal" polarizers (a pplicable for ACS, b ut not for 
NICMOS) is discussed bv lCapetti eFafl (l2fJ0l . 

The instrumental counts per second measured through 
each polarizing element are used to define an observed 
intensity vector of the form a = [Ji, 12,13]. The Stokes 
parameters, which also defined a vector b — [I,Q,U], 
are used to calculate p and 9f (see Eqs. [5] and [3]). The 
two vectors a and b are simply related to each other by 
the linear expression [C]b = a, where [C] is a matrix de- 
scribing the characteristics of the k th polarizer, namely 
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Fig. 6. — The results using the previous calibration. The thick solid lines show the ground based results and the grey shaded areas show 
the uncertainties associated with those ground-based results. To avoid cluttering the plot, error bars are only shown for the cases where 
they do not overlap the ground based results. The vertical solid lines mark the radii of the 1 st , 2 nd and 3 rd dark minima of the PSF. The 
values of 9 are in the celestial reference frame (position angle east of north). The values of 9 for HD 331891 bear no meaning and have 
therefore been excluded. 



TABLE 4 

The NIC2 Calibration from Multiple Orientations 



Target 


Polarizer 


h 


h 


h 


p% 




VR 84c 


POL0L 


22357.5 


22855.3 


22578.1 


1.7 ±0.3 


115 ± 2 




POL120L 


22914.5 


22751.2 


22721.1 


0.8 ±0.3 


126 ± 6 




POL240L 


22775.4 


22650.3 


23156.0 


1.6 ±0.5 


122 ± 6 


VSS VIII-13 


POL0L 


34362.8 


34586.6 


34804.0 


1.1 ±0.3 


79 ± 10 




POL120L 


34733.7 


34816.8 


34635.1 


0.5 ±0.2 


85 ±22 




POL240L 


34969.2 


34487.8 


35190.2 


1.3 ±0.3 


105 ± 6 


HD 331891 


POL0L 


33817.0 


33881.5 


33636.2 


0.6 ±0.6 






POL120L 


34154.7 


34081.5 


33929.6 


0.7 ±0.2 






POL240L 


34244.6 


34063.8 


34054.5 


0.4 ±0.2 





Note. — All results have been extracted from the apertures with a radius 
of 7.6 pixels (0"58) and combines the data from the individual polarizers at the 
separate roll angles. Ii , I2 and 73 are in units of counts per second, averaged 
across the dithers, per polarizer, per orbit. 



the actual orientation (in radians) of the polarizer (4>k), 
the fraction of the light transmitted in the parallel di- 
rection (tk), and the fraction of light transmitted in the 
perpendicular direction (the "leak", The polarizer 
efficiency (77^-) is given by (1 — 'fe)/(l + h)- The linear 
expression [C]b — a can be solved for (I, Q, U) using LU 
decomposition. We can then determine p and 9f using 
Equations [2] and [3l 

P = 100%x^7^ (2) 

6>f = -arctan— (3) 
1 2 Q w 

In Equation [3] a 360° arctangent function is assumed. 
In addition, the orientation of the frame has to be sub- 
tracted from 6 f in order to retrieve the celestial position 



angle 9. This process has been coded into the IDL rou- 
tine polarize. pro 10 . This routine takes the three polarized 
images, calculates (I,Q,U), and then produces two di- 
mensional maps of p and 8 f . The la dispersions of the 
count rates around the mean (used to define S/N CT ) have 
also been used to determine the uncertainties in p and 9. 
As shown by SA99, the variance and covariance in the 
(Q,U) plane define the uncertainties in p (o~(p)) and 9 
(a(9)). 

In order to demonstrate the limiting accuracy of the 
previous (Cycle 11) polarimetric calibration, we have ap- 
plied the Table [3] coefficients to our newly acquired po- 
larimetric standard star observations. Figure [5] demon- 
strates that only at a few roll angles does this reproduce 
the ground-based results. In all cases, however, we can 

10 http:/ /www. stsci.edu/hst/nicmos/tools/polarize_tools. html 
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Fig. 7. — The results using the new calibration. The same convention is used as for Figure [6] All results are now consistent with the 
ground-based data. The values of 9 for HD 331891 bear no meaning and have therefore been excluded. 



TABLE 5 
The NIC2 Instrumental 



Target 


Polarizer 




Qnic 


Unic 


Qstd 


U s td 


Pins% 


@ins ( ) 


VR 84c 


POL0L 


44574.1 


-490.7 


-573.8 


-172.2 


-529.9 


0.7 


94 




POL120L 


44315.3 


-112.7 


-355.4 


-171.2 


-526.8 


0.4 


126 




POL240L 


44294.8 


-309.4 


-613.7 


-171.1 


-526.6 


0.4 


106 


VSS VIII-13 


POL0L 


68183.2 


-680.6 


281.6 


-566.8 


-252.4 


0.8 


141 




POL120L 


67523.4 


-325.9 


52.6 


-561.3 


-249.9 


0.6 


116 




POL240L 


67563.3 


-787.7 


-439.9 


-561.7 


-250.1 


0.4 


110 


HD 331891 


POL0L 


66628.9 


85.1 


-399.8 


0.0 


0.0 


0.6 






POL120L 


66209.9 


-315.8 


297.2 


0.0 


0.0 


0.7 






POL240L 


66142.9 


249.9 


-121.8 


0.0 


0.0 


0.1 





Note. — The results from the NIC2 instrumental polarization tests. All results have 
been extracted from the apertures with a radius of 7.6 pixels. 



see that the results are stable outside of the 2 nd Airy 
maximum. 

4.1. The Calibration 

It is now that we can begin to exploit the power of ob- 
serving each standard at multiple spacecraft roll angles; 
the uncertainties due to the relative hlter transmissions 
can be removed as the same polarizers have been used at 
three distinct orientations, i.e., each polarizer has sepa- 
rately, and independently, collected all of the necessary 
polarimetric data. By using polarimetric data gathered 
in this fashion, the method for calculating p and 9 can use 
invariant transmission coefficients for a given polarizer. 
In this case the resulting polarization is independent of 
t k and is only affected by 4> k (plus the ORIENTAT) and 
Ik- Due to observation scheduling constraints imposed 
by HSTs 2-Gyro operating mode (in place during Cy- 
cle 15) we were unable to obtain our calibration data at 
ideal differential celestial position angles (see Table [T] for 
as-executed orientations). Hence, though obtained with 
three roll angles, we cannot analyze these data as one 



would through a single perfect polarizer at ideal orienta- 
tions; the method described in §[4] must still be followed 
in order to derive the Stokes parameters. 

The quantitative results obtained from polarimetric 
analysis of the apertures enclosing the 2 nd Airy maxima, 
from the three roll orientations using the Cycle 11 derived 
calibration coefficients (Table |3|) , are shown in Table |H 
These results do not precisely recover the ground-based 
results listed in Table [2j With the unknowns (tk) nulli- 
fied by the three orientations, we assume these residuals 
to be an instrumental polarization (pi ns )- 

Since (Q, U)i ns is in the instrumental frame, the de- 
rived pi ns is fixed to the co-rotating spacecraft frame. 
Therefore, simply subtracting the observed residual po- 
larizations (found in the unpolarized standard) from the 
polarized standards is not appropriate; pi ns must be de- 
termined in Stokes space. For this, as we are calibrating 
NICMOS against the ground-based results, we determine 
(I, Q, U) (and the associated uncertainties) for the stan- 
dard stars, (/, Q, U) s td- However, / is dependent on the 
nature of each individual observation and the intrinsic 
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TABLE 6 

Adjusting the Transmission Coefficients. 



Target 


PA 


*o 


*120 


P (%) 


e 


VR84c 


-178.6° 


0.8827 


0.8346 


1.2 ±0.4 


126 ±3 




-61.4° 


0.8804 


0.8369 


1.3 ±0.4 


126 ±5 




58.4° 


0.8883 


0.8447 


1.3 ±0.4 


127 ±7 


(B06) 


-50.8° 


0.8914 


0.8372 


1.2 ±0.1 


125 ± 7 


(B06) 


-160.8° 


0.8897 


0.8397 


1.2 ±0.1 


124 ± 6 


VSS VIII-13 


179.6° 


0.8737 


0.8335 


0.9 ±0.4 


103 ± 12 




113.6° 


0.8771 


0.8369 


0.9 ±0.4 


103 ± 13 




53.6° 


0.8916 


0.8425 


0.8 ±0.3 


102 ± 8 


HD 331891 


-86.9° 


0.8804 


0.8391 


0.3 ±0.4 






-22.4° 


0.8860 


0.8302 


0.5 ±0.5 






37.5° 


0.8849 


0.8391 


0.4 ±0.4 




(B06) 


-109.6° 


0.8629 


0.8269 


0.3 ±0.1 




(B06) 


0.6° 


0.8718 


0.8334 


0.3 ±0.1 





Note. — Determining the parallel transmission coefficients 
that reproduce the ground-based results (taken from a 7.6 pixel 
aperture). 



polarization of the target, so we must calculate (Q, U) s td 
as a function of the / observed by NIC2 (I n ic)- It then 
follows that we can determine (Q, U)i ns from Equation^ 
(where all vectors have been rotated into the celestial 
plane). We can then use Equations [2] and [3] to deter- 
mine (p,9)i ns . The results of this are presented in Ta- 
ble [51 which shows the raw Stokes parameters derived 
from following this method. The instrumental polariza- 
tion derived from an unweighted average of our observa- 
tions is 0.6 ± 0.1% at 116 ± 15°. 

(Q,U) ins = (Q,U) nic -(Q,U) std (4) 

Attention can now be turned back to the data from the 
single roll angles. With the instrumental polarization (in 
Stokes space) determined (and subtracted), the only pa- 
rameter affecting the deviation from the ground-based 
results is tk- The values of to and £120 that reproduce 
the ground-based results can then be determined numer- 
ically. As in the previous calibrations, £240 is held at 
0.9667 as it is the long polarizer with the highest effi- 
ciency. The other values of tk can be determined with 
respects to £240- 

As this process was performed by B06 on archival data, 
we can also include the values of tk derived by those 
authors. We limit these additional data to the polari- 
metric standards observed in the post-NCS era. These 
data comprise of two epochs for VR 84c (identified as 
CHA-DC-F7 by B06) and two epochs of data for HD 
331891. Table [5] presents the resulting values of to and 
ti2Q including the derived values of p and 9. The un- 
weighted average values of t and ti 2 o are 0.882 ± 0.008 
and 0.837 ± 0.004 respectively. 

We apply these newly re-derived coefficients to ascer- 
tain (p,9) for our target stars with the results shown 
in Figure [7] (comparable directly to Figure [6]) and find 
the new calibration produces results consistent with the 
ground-based determinations in all cases. As an addi- 
tional check we also apply the updated transmission co- 
efficients to the data presented by B06 that suggested the 
original instrumental polarization. Figure [5] plots the re- 
sults from these re-analyzed data (it is similar to Figure 
2 from B06). As can be seen, re-analysis with the now 
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Fig. 8. — Checking the new calibration against the data that 
originally suggested an instrumental polarization. The errors bring 
all consistent with the ground-based results. The values of 9 for 
HD331891 have been omitted for clarity, and because they bear no 
meaning. 

improved determinations of the polarizer transmissivitics 
produce (p, 9) results that are also now consistent with 
the ground-based data. The typical uncertainties asso- 
ciated with the profiles in Figure [5] are a p « 1% and 
<jg 10°. The Cycle 11 data were not gathered using 
the additional dither points employed by this project; the 
polarimetric analysis performed on the archival data was 
improved through our updated calibration. 

5. DISCUSSIONS 

We have used our data to perform a detailed investi- 
gation of the NIC2 instrumental polarization. This work 
was motivated, in part, by recent attempts to observe 
objects that have intrinsic polarizations of less than 5%. 
The previous calibrations were of sufficient accuracy to 
complete studies of more highly polarized targets but 
could not be used for high accuracy polarimetry of tar- 
gets with low intrinsic polarization. By obtaining high 
S/N observations of one unpolarized and two polarized 
standard stars, at three well separated roll angles, we 
have been able to determine the level and orientation of 
the NIC2 instrumental polarization. By subtracting the 
instrumental polarization, in Stokes space, and improv- 
ing upon the previous determination of the parallel trans- 
mission coefficients, we have successfully reproduced the 
polarimetric parameters of three standard stars observed 
from the ground. These results were consistent for each 
standard observed at three different instrumental orien- 
tations in the celestial frame. Our results make this the 
most comprehensive polarimetric calibration of NIC2 to 
date (and quite likely the last) for direct imaging. We 
suggest this quantitative characterization of the NIC2 
polarimetric system supersede the previous post-NCS era 
polarimetric calibrations. 

As the techniques used here require extremely pre- 
cise photometry, there are many potential sources con- 
tributing to the uncertainties in determining p and 9. 
Highly precise photometric repeatability is limited by 
temporal, thermal, and metrologic optical instabilities 
in the J?5T'+NICMOS+ instrumental system affecting 
the PSF, as well as spatial quantization and perfor- 
mance defects in the NIC2 detector and readout elec- 
tronics. The measurement dispersions introduced by 
these (and perhaps other systematic) effects are car- 
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lied through the SA99 analysis to determine <j(p) and 
cr{9). However, as we have been able to determine re- 
vised values for to and £120 in 13 cases, we can assign 
uncertainties to these coefficients. By propagating these 
uncertainties through the polarimetric analysis we find 
that they dominate over the errors derived following 
SA99. For example, applying the new coefficients to 
VR84c we find p = 1.3 ± 0.4,0 = 126 ± 5° following 
SA99, but by propagating the uncertainties in tk we find 
p = 1.3(+0.7, -0.6)%, 9 = 126(+10, -24)°. 

The remaining uncertainties in the derived transmis- 
sion coefficients do indeed dominate the statistical un- 
certainties from the photometry. What affect do these 
uncertainties have on the observationally determined val- 
ues of pi ns and 0i ns ? Could these uncertainties be re- 
sponsible for the observed instrumental polarization? As 
can be seen from Table 21 deviations between I\ , 1% and 
-Z3 (that in an ideal detector would be invariant) are no 
greater than ~ 2% which corresponds to a la error in 
Pi ns of only 0.1%. However, a more accurate way to 
determine pi ns requires observations of an unpolarized 
standard star across the whole detector area to map any 
field dependence in the polarization. As our results have 
been derived from dithered observations, if there is a field 
dependence in the polarization, then that could mimic an 
instrumental polarization on our data. Such a study can- 
not be done with the current data, so it is not possible to 
assert pi ns at such a level. Instead, it is more appropri- 
ate to present the 0.6% value as an instrumental upper 
limit; pi ns could indeed be zero. 

What then happens to the rest of our analysis if pi ns 
is indeed zero? To test this, we have re-derived the val- 
ues of tk without subtracting pi ns in the (Q, U) plane. 
Taking this approach we find to — 0.884 ± 0.008 and 
t%20 = 0.838 ± 0.003 which is entirely consistent with the 
results taken from the instrumentally subtracted (Q, U). 
So, again we find that the errors in tk dominate the un- 
certainties in the measured values of p and 9. As these 
uncertainties have been derived from 13 separate obser- 
vations, that are all subjected to the photometric uncer- 
tainties themselves, we are left with the understanding 
that this calibration is limited by the photometric re- 
peatability of individual observations. This includes all 
of the known affects presented in § [2] 

When comparing the new parallel transmission coeffi- 
cients with the old it can be seen that t\2o is consistent 
with both calibrations; only to has needed to be changed, 
albeit "only" by ~ 0.5%. Applying this correction to 
previous studies of highly polarized targets will therefore 
not alter the results by any significant amounts. How- 
ever, when applied to targets with low polarization this 
"tweak" has a tremendous affect. 

5.1. Recommended Observing Strategy 

In the cases where the target is spatially resolved, and 
does not contain a strong point source, then observers 
should bin their data to 3 x 3 pixels (in order to meet 
the diffraction limit of the instrument). Dithering can 
be used in the usual way to increase spatial resolution 
as well as to mitigate bad pixel affects. However, signal 
to noise will still limit the accuracy of p and 9 determi- 
nations, especially in targets whose surface brightnesses 
vary significantly on small spatial scales, i.e., comparable 
to, or smaller than, the diffraction limit. 



TABLE 7 
The New Polarimetric coefficients 



Polarizer 


4>k (°) 


Vk 


k- 


tk 


POL0L 


8.84 


0.7313 


0.1552 


0.882 ± 0.008 


POL120L 


131.42 


0.6288 


0.2279 


0.837 ± 0.004 


POL240L 


248.18 


0.8738 


0.0673 


0.9667 



Note. — The new NIC2 polarization coefficients as 
derived by this program. 



The case of a spatially unresolved source is not so 
straightforward. As explained above, the uncertainties in 
tk do dominate the accuracy of p and 9. However, as this 
study has shown, non-repeatabilities in aperture pho- 
tometry (in the absence of well-dithered observations) 
suggest that unidentified outliers can readily bias polari- 
metric analysis independent of the intrinsic accuracy of 
the polarimetric calibration. The data have been care- 
fully inspected to identify photometrically deficient pix- 
els by comparing each individual profile to the dispersion 
in the photometry measured from ~ 12 pointings. This 
study has shown that ~ 10% of all target pointings (even 
those avoiding known photometrically deficient areas of 
the detector like the coronagraphic hole) will include un- 
correctable bad pixels that affect the measured polariza- 
tion. A method to robustly detect, and reject, photo- 
metric measures degraded by such effects is to perform 
a target raster with a dither pattern with three or more 
pointings. The step size should be more than three times 
greater than the FWHM of the PSF to avoid persistence. 
Encircled energy (intensity) profiles should be derived in- 
dependently from each observation (dither point). With 
a sufficiently large number of dither points, a sigma clip 
can then be used to determine whether an individual pro- 
file is deviant or not. The greater number of dithers used, 
the more accurate this clip will be. For small numbers 
of dither points (not optimally recommended) medianing 
can be used to coarsely reject outliers. The average of 
the remaining profiles can then be used to derive p and 
9 from the coefficients presented in Table [7J 

This study has shown that the affects of the IPRFs are 
largely mitigated by effective dithering. However, the un- 
certainty in the calibration of tk ultimately dominates the 
accuracy of (p, 9) polarimetric analysis for sources with 
intrinsically low polarization, if observed with a sufficient 
number of optimally placed dither positions; observers 
need not be concerned by the IPRF. Since observations 
through the three polarizers are not simultaneously ob- 
tained, temporal instabilities in the fiST+NICMOS PSF 
are still a concern and the effects are readily seen in Fig- 
ures [2 [6] and [Jj Outside of an aperture of radius 0758, 
the PSF affects are seen to be alleviated. Inside this ra- 
dius, the errors in p and 9 rise rapidly and it will be left 
to the discretion of the observer to weigh their required 
accuracy to the possible results of beam-depolarization. 

6. CONCLUSIONS 

In a non-ideal imaging polarimeter, such as NICMOS, 
it is essential to observe several polarimetric standards, 
at three well separated position angles through each po- 
larizer, in order to full characterize the instrumental po- 
larization. The additional roll angles remove the un- 
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Fig. Al. — The NIC1 radial profiles. For each target, the encircled count rate for each polarizer and celestial orientation, is shown. 



certainties in the (unknown) parallel transmission coeffi- 
cients and allow rotation of any instrumental polarization 
with respects to the equatorial frame. Using this tech- 
nique we have placed an upper limit to the NIC2 instru- 
mental polarization of 0.6%. With a known value for the 
Stokes / parameter, the instrumental polarization can be 
transformed into the Q, U plane and be subtracted. New 
parallel transmission coefficients can then be determined 
numerically by comparing p and 9 with that of the apri- 
ori well-determined calibration standards. Following this 
approach we have determined the to and ti20 coefficients 
to be 0.883 ± 0.004 and 0.837 ± 0.004 respectively. As 
with the previous calibrations of the NIC2 polarimeter, 
we held t24o constant at 0.9667. The t\2o coefficient is 
consistent with the previous calibration, but our knowl- 
edge of t has been improved, resulting in a change in 
its previously determined value by ~ 0.5%. Such a small 
change in to does not warrant the re-analysis of previ- 
ous NIC2 imaging polarimetry data of highly polarized 
(p > 5%) targets, but is significant for targets with in- 
trinsically low polarization fractions. 

As we use 13 determinations of tfc, we are able to assign 
la uncertainties to the calibration coefficients. Propa- 



gating these uncertainties through the polarimetric anal- 
ysis, we find that they dominate all other sources of er- 
ror. Applying these adjusted values of tk to this (and 
the archived) calibration data, we find that NIC2 is now 
capable of confidently detecting polarizations at a level 
of 1.0%. The uncertainties associated with such mea- 
surements are ±0.6% and ±15° inp and 6, for sources at 
the w 1% level of intrinsic polarization. This is the first 
time that such a level of accuracy has been achieved with 
the NICMOS polarization calibration. This improved 
calibration opens a new domain for observational inves- 
tigations with HST by enabling very high precision po- 
larimetry of intrinsically very low polarization sources. 
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APPENDIX 
THE NIC1 CALIBRATION 

For the NIC1 observations, the individual pointings for each polarizer were not dithered around the 12 point pattern 
used for the NIC2 observations. However, the observations from each separate polarizer were dithered so as to avoid 
persistence. The pointings were also placed on an area of NIC1 where the quantum efficiency (QE) gradients are 
particularly shallow and where the QE is reasonably high. This area was also sufficiently far from any "grotty" pixels 
and sufficiently far from the edge of the detector. 

As there are only single pointings per polarizer, the data analysis for the NIC1 data is less complicated than for the 
NIC2 case. The extracted radial profiles from each polarizer (and pointing) at each roll angle are shown in Figure IAT1 
For these profiles we can only define the S/N using the photon noise. As there are only single pointings in each NIC1 
polarizer, it is impossible to perform the iterative 2a clipping technique employed for NIC2. Therefore all NIC1 profiles 
must be included in the analysis. Although there is no statistical way to tell if the profiles in Figure [ATI are robust, 
we can at least see that all profiles are consistent with each other. This is likely due to the care taken in avoiding 
uncorrectable bad pixels. 

At the radius of the 3 rd NIC1 dark Airy minimum (given by an aperture with r = 7.0 pixels, assuming a central 
wavelength of 1.05/um and a 2.281m aperture) the respective signal to noise ratios for VR84c, VSS VIII-13 and 
HD331891 are 257, 180 and 330 when averaged across the three polarizers. The theoretical NIC1 achievable accuracies 
are approxim ately 1.5% wit h an u ncertainty of ±0.2% and 1.0 ± 0.3%. 

At 1.21jum. lWhittet et al.l(|1992t ) give VR84c a polarization of 3.19±0.05% at 118±1° and VSS VIII-13 a polarization 
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TABLE Al 

The NIC1 Calibration from Multiple Orientations 



Target 


Polarizer 


h 


h 


h 


p% 




VR84c 


POLOS 


64292.5 


68734.3 


65147.3 


4.1 ± 0.1 


116.5 ±0.5 




POL120S 


67785.9 


64160.6 


66039.8 


6.47 ±0.06 


131.6 ±0.6 




POL240S 


66771.2 


64206.3 


68756.5 


5.1 ± 0.2 


122.8 ±0.2 


VSS VIII-13 


POLOS 


31156.2 


32770.7 


32472.4 


3.1 ± 0.2 


91 ±1 




POL120S 


32393.0 


32412.3 


31762.6 


3.0 ±0.1 


83 ± 2 




POL240S 


33123.9 


31812.0 


32638.4 


2.9 ±0.2 


94 ±2 


HD 331891 


POLOS 


107931 


107057 


109394 


1.3 ± 0.1 






POL120S 


108864 


109376 


108867 


0.63 ± 0.04 






POL240S 


108221 


108756 


109656 


1.06 ± 0.06 





Note. — All results have been extracted from the apertures with a radius of 
7.0 pixels and combines the data from the individual polarizers at the separate roll 
angles. 



TABLE A2 

The New NIC1 Polarimetric coefficients. 



Polarizer 


4>k (°) 


Vk 


h 


tk 


POLOS 


1.42 


0.9717 


0.0144 


0.7760 


POL120S 


116.30 


0.4771 


0.3540 


0.57 ±0.01 


POL240S 


258.72 


0.7682 


0.1311 


0.69 ±0.01 



Note. — The new polarization coefficients, as 
derived by this program, for NIC1 polarimetry. 



of 2.21 ±0.3% at 96±1°. As there is no data for HD331891 at this wavelength we assume that this target is unpolarized. 
The polarizations and uncertainties are carried through the Serkowski correction and are presented in Table [5J 

Except for the 2a clipping, the NIC1 data was treated exactly like the NIC2 data. The calibration data from 
the multiple orientations are presented in Table [ATI The NIC1 instrumental polarization was then determined to be 
1.3±0.9% at an orientation of 151 ±18°. This instrumental polarization was subtracted in (Q, U) space and the parallel 
transmission coefficients were re-derived on an orbit by orbit basis. The coefficients, and uncertainties, that bring all 
the data consistent with the ground-based data are t\2o = 0.57 ± 0.01 and £240 = 0.69 ± 0.01. All the coefficients to 
be used in NIC1 polarimetric observations are shown in Table lA"2l Applying these new calibration coefficients to the 
original data brings all polarimetric standards consistent with the ground-based results. 

Compared with the previous NIC1 calibration coefficients (£120 = 0.5934 and £240 = 0.7173) the coefficients derived 
here are lower. In the previous instance the coefficients were derived from observations of an unpolarized standard, and 
not from data at three well separated roll angles. However, the data used in this calibration have not been dithered over 
the detector (as in the case of the NIC2 data); instrumental affects will still be present. The fact that the uncertainties 
in the newly derived coefficients for NIC1 are larger than for the NIC2 case, reflects this finding. Propagating the 
coefficient uncertainties through the polarimetric analysis results in polarization measures of 5(+6,-2)% at 123(+24,- 
12)°for VR84C, 3(+6,-2)% at 94(+85,-27)°for VSS VIII-13 and 3(+5,-3)% for HD331891. We therefore recommend 
that NIC1 only be considered for observations of targets that are postulated to have an intrinsic polarization of greater 
than 4% per resolution element. The uncertainties in polarization measures from such a target will, however, be large. 
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